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STABILITY ESTIMATES FOR THE INVERSE BOUNDARY 
VALUE PROBLEM FOR THE BIHARMONIC OPERATOR 
WITH BOUNDED POTENTIALS 


ANUPAM PAL CHOUDHURY AND VENKATESWARAN P. KRISHNAN 

Abstract. In this article, stability estimates are given for the deter¬ 
mination of the zeroth-order bounded perturbations of the biharmonic 
operator when the boundary Neumann measurements are made on the 
whole boundary and on slightly more than half the boundary, respec¬ 
tively. For the case of measurements on the whole boundary, the stabil¬ 
ity estimates are of ln-type and for the case of measurements on slightly 
more than half of the boundary, we derive estimates that are of In ln- 
type. 


1. Introduction 

Let 17 C M n , n > 3 be a bounded domain with C°° boundary and consider 
the following equation: 

B q u := (A 2 + q)u = 0 in 17, q € L°°(17). 

Let the domain of B q be 

D(B q ) := {u € H\n) : u\ an = Au\ m = 0}. 

We consider the following space for the potential q: 

(1) Qm '■= {q ■ supp(g) C 17, and ||r?||(f 2 ) A M for some M > 0}. 

We will assume that for all q € Qm, 0 is not an eigenvalue for B q on the 
domain D{B q ). Then given (/, g) € H 7 / 2 {dVt) x R 3//2 (<9I7), there is a unique 
solution to the boundary value problem: 

(2) B q u = 0, u\an = f, Au|an = g- 

The boundary conditions are called Navier conditions [6] and we define the 
Dirichlet-to-Neumann map J\f q for this operator by 

M q :R 7/2 (<917) x R 3 / 2 («917) -»• R 5 / 2 («917) x H 1 / 2 ^ 7) 

., . ,du, d(Au). . 

(3) (f,9)->(-Q^\dSh—Q^—\dn)> 

where u € 7) is the unique solution to ([2]). 

We are interested in the inverse problem of determining q from J\f q . The 
uniqueness question of determination of q from J\f q was answered in mm 
and recently in rnmm where they showed that unique determination 
of both zeroth- and first-order perturbations of the birharmonic operator is 

l 
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possible from boundary Neumann data. We note that the papers DUE! 
also show unique determination of the first order perturbation terms from 
Neumann data measured on possibly small subsets of the boundary. 

In this paper, we consider the stability question for the determination of 
q from J\f q for the operator B q . That is, whether one can estimate perturba¬ 
tions of q from perturbations of the Neumann data J\f q . To the best of the 
authors’ knowledge, stability estimates for inverse problems involving the 
biharmonic equation has not been obtained earlier, and the purpose of this 
paper is to investigate it. We prove a stability estimate of ln-type for the 
case when the Neumann data is measured on the whole boundary. We then 
prove a stability estimate of lnln-type when the Neumann data is measured 
on a part of the boundary that is slightly more than half the boundary. 

Our strategy for proving stability estimates follows the methods intro¬ 
duced by Alessandrini in [T] using complex geometric optics (CGO) solu¬ 
tions where a ln-type stability estimate is proved for the Calderon inverse 
problem [3], and by Heck-Wang in [8] where a lnln-type stability estimate 
is proved for the Calderon inverse problem when the Neumann data is mea¬ 
sured on slightly more than half of the boundary. CGO solutions were in¬ 
troduced by Sylvester and Uhlmann in the fundamental paper HU to prove 
global uniqueness for the Calderon inverse problem. The method in Heck 
and Wang combines CGO solutions and techniques of [2] with an analytic 
continuation result of Vessella |16j . Stability estimates for several inverse 
problems have been obtained in recent years. Apart from the works in Ei 
already mentioned, we refer the reader to D3G1I3Q] for stability estimates 
involving the Calderon inverse problem and inverse problems involving the 
Schrodinger or magnetic Schrodinger equation. 

2. Statements of the main results 

We now state the main results of this paper. We first consider stability 
estimates for full boundary measurements and then prove stability estimates 
when only partial boundary measurements are available. 

2.1. Results for full boundary measurements. Consider the following 
norm on H a (dQ) x (for simplicity we will denote this space by 

H a ’P(dn)): 

(4) \\(f,g)\\H°‘’P(dn) = 11 /1 In** (an) + llfl , lln/ 3 (<9n) f° r (/>S0 € H a - 6 (dLl). 

Define: 

Ml = s u p{||A / -,(/, 5 )||^,i (an ) : \\(f,9)\\ H H (da) = 1} 

where M q (f,g) is defined in ©. 

Theorem 2.1. Let D C M n ,n >3 be a bounded domain with smooth bound¬ 
ary. Consider Equation (J2J) for two potentials q\,q 2 € Q■ Let M qi and J\f q2 
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be the corresponding Dirichlet-to-Neumann maps measured on dfl. Then 
there exists a constant C = C(fl,n, M) such that 

hi ~ Q2\\ 2 H -i {n) <C [\\M qi - M q2 \\ + | InHA/'g, . 

2.2. Results for partial boundary measurements. Now we consider 
the problem of estimating perturbations of q, when the Neumann data J\f q is 
measured on a subset of dfl that is slightly more than half of the boundary. 

Before stating the result, we introduce the following notation. Let a G 
S 71 ” 1 be a unit vector and e > 0 be given. Let v{x) denote the outer unit 
normal at x G dfl. We define 


(5) = {x G dfl, a ■ v(x) > e}, dfl- >£ = dfl \ dfl 

(6) dfl + = {x € dfl, a ■ v{x) > 0}, = dfl \ 


Now the partial Dirichlet-to-Neumann map is defined as 


M q :iL2’2(<9fl) -A- #2’2 (<9D_ i£ ) 

(/)<?) ->• {dvu\ a n_ e ,d„(Au)\ a n_ e ), 


where u € H i {fl) is the unique solution to ([2]). As before, we define the 


norm of Af q as 


Wq\\= SUp{||A/'g(/,, 




= 1 } 


We have the following stability estimate with partial boundary measure¬ 
ments. 


Theorem 2.2. Let D C M n ,n >3 be a bounded domain with smooth bound¬ 
ary. Consider Equation m for two potentials qi,q 2 € Q. Let J\f qi and N q2 
be the corresponding Dirichlet-to-Neumann maps measured on dfl- t£ . Then 
there exist constants C = C(fl,n,M,e), K and 9 > 0 such that 


19i -®||H-i(n) < ll-Mfi -AT q2 \\ + ( — In | lnHAAg 


J q i 



o 
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3. Preliminary results 

We use the following result from [IIlEI- 

Proposition 3.1. [11, Prop. 2.2] (Interior Carleman estimates) Let q G 
Qm and p = x ■ a, |a| = 1. There exists an 0 < ho = ho(n,M) <C 1 and 
C = C(n, M) > 0, where n is the dimension and M is the constant in (jT]) 
such that for all 0 < h < ho <C 1 and u € Cf°(fl), we have the following 
interior estimate: 

\\ef' h h i B,e-^ h u\\ L , {n) > L||„|| ff . i(n) , 

This result is based on a Carleman estimate proven in 031 - 
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Proposition 3.2. [Ill Prop. 3.2] (Boundary Carleman estimates) Let q € 
Qm and ip = x ■ a, |o| = 1. Let be as in (JB]). There exists an 0 < ho = 
ho(n,M) <C 1 and C = C(n,M ) > 0, where n is the dimension and M is 
the constant in <tu> such that for all 0 < h < ho <C 1 and u G D(B q ), we 
have the following estimate involving boundary terms: 

||e _¥,//l /i 4 5 gu|| L 2 (n) + h 3/ ~\\s/-a ■ ve~ v/h d u {-h 2 Att)|| L 2 (90 _) 

(7) 

+ h^V^^e-«/ h d„u\\ L2{dn _) > ±(h 2 \\e-*/ h u\\ HUa) 

+ hW\Wa ■ ve h 2 Au)\\ L 2 ^ + ) + h°'~\\^a ■ ve v ^ h d v u\\L 2 ^Q n+ ^j. 

Using estimate of Proposition l3.ll the following result is proven in [11. 12] 
which we will require in what follows. 

Proposition 3.3. [Ill Prop. 2.4] There exists an ho = ho(n,M) > 0 and 
C = C(n, M ) > 0, where n is the dimension and M is the constant in (JXj) 
such that for all 0 < h < ho <C 1, there exist solutions u(x, £; h) G H l (Ll) to 
B q u = 0 in Ll of the form 

u(x, £; h ) = e~~ (1 + hr(x, (; h)), 

with £ G C n satisfying £ ■ £ = 0, |i2e(£)| = |/m(£)| = 1 and |]r||#4 < Ch 2 . 

scl (£2) 

We note that the estimates on ho and r are independent of the potential 
q G Qm- 

For proving stability estimates with partial data, we require the following 
result due to Vessella [16] . 

Theorem 3.4. [16], Theorem 1] Let Q C M n be a bounded open connected 
set such that for a positive number ro the set £l r = {x G Ll : d(x, dQ) > r} 
is connected for every r G [0, ro]. Let E C be an open set such that 
d(E,dQ ) > do > 0. Let f be an analytic function on with the property 
that 

\D a f(x)\ < ^ for x G n, a € (N U {0}) n , 
where A, C are positive numbers. Then 

/ \7i(|E|/M) 

|/(x)|<( 2 U) 1 -^™ n ^sup|/(x)|J 

where |i£| and |ST| denote the Lebesgue measure of E and Ll respectively, 
7 i G (0,1) and 71 depends only on do, diam(il), n, ro, A and d(x,dQ). 
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(8) 


We also require the following Green formula: 

J(B q u)vdx — J uB*vdx = j d u (Au)v dS + j d u u(Av ) 


dS 


an 


an 


— J (A u)d u v dS — J u(d v {Av)) dS 1 . 
an an 

4 . Stability estimates with full boundary measurements 
In this section, we prove Theorem 12.11 

Proof of Theorem \2.1[ We start with the Green formula © and let q = qi 
and u = u\ — U 2 and v € H' [ (Q) is such that B* n v = 0 in Q. Here u\ and u 2 
are solutions to (J2J) for q replaced by q\ and q 2 - Then we have 

(9) j (92 — qi) u 2 vdx = J d v (A(ui - u 2 ))v dS + J d v {u\ - u 2 )( Av) dS 1 . 


an 


an 


Using Proposition 13.31 we have solutions to B q2 U 2 = 0 and B* qi v = 0 of the 
form 


( 10 ) 

( 11 ) 

where 


■■<1 


v(x,(i',h)=e h (1 + hr^x,Cl,h)), 


<2 


u 2 {x,( 2 -,h) = e h (l + hr 2 (x,C 2 ',h)), 


Cl = f + V i -'> 2 x ,3+i “’ 


b=- h 4 + fi-P^- m . 

with a and (5 are unit vectors in M n with a, /3 and £ are mutually perpendic¬ 
ular vectors and h is such that h < ho and 1 —is positive. Substituting 
U 2 and v into the left hand side of we get, 

(12) _ 

I (92 ~ qi)u 2 vdx = (q 2 - 9 i)(£) + / (92 - qi) e~ lx< (hrJ + hr 2 + h 2 nr 2 ) dx. 
J Jsi 

n 

Calling the second term on the right hand side of the above equation as /, 
we have the following estimate. 

\I\< / 1 92 - qi\{h\n\ + h\r 2 \ + /i 2 |ri||r 2 |) dx 
Jn 

< C(h\\ri\\ L 2^ + /i||r 2 || L 2(Q) + /i 2 ||ri|| L 2(Q)||r' 2 || L 2(Q)) 

< Ch since /i< 1. 


( 13 ) 
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Now consider the right hand side of (j9]). We have 

| [ d u (A(ui — u 2 ))vdS + f d u (ui - u 2 )(Av)dS\ 

Jan Jan 


< 


i an 

f \d u (A(m - u 2 ))v\ dS + f \d u (u\ - u 2 )(Av)\dS 
Jan Jan 




< \\d u (A(ui — ^2 ) ) 11Z, 2 (SO) 11 ^ 11Z. 2 (0O) + \\d u {ui - ^2)||L2(dO)l|At’||L2( 9 Q) 

< C(\\d u (A(u 1 - w 2 ))Ilz, 2 (ar2)ll^llizi(si) + _ M 2)||L 2 (an)l|Au||#i(n)) 

< C(\\d u (A(ui — u 2 ))\\ L 2( d n) + \\d u (ui — rt2)||i2(an)) + ll^ v ll_ff 1 (n)) 

< C7(||^ / (A(«i -u 2 )),d u ( Ul -«2)ll^,§ (an) )(l|u|lHi ( n) + H A «lliP(n))- 

which again is 

= cm qi -Af q 2 )(f, i 

< C\\N qi -Nq 2 \\\\Ui 

< C\\M qi - Afq 2 \\{\\ u 2 \\m(ci) + \\Au 2 \\ H 2 (n)) (IMIfD(fi) + ll^Olff 1 ^))- 

We have the following estimates for ||v||#i(n) an d ||Au||^i(q). In these esti- 

ix-Cj 2 R 

mates, we use that fl C B(0,R) for R > 0 fixed. Then |e ^ | < e~, since 
101 = 2 for j = 1,2. 

n 

i3>Cl ^—> ix-Ci 1 

IMIifi(n) - H e ~( 1 + hr i)\\L\n)+2^W he ~d Xk ri + -Cifee~(l + hn)||i,2 (n) 


5 1 


fc=i 


+ 


n / , 

<e~(l + h\\ri ||//4 i ( f2 )) + ^2 ( 2e ~ ( ll ri 

k =1 ' ' 

2 R . r\ Cy 2 R . r\ Cy 2R 

< Ce h (1 + h 2 ) + — e h (1 + h 2 ) < — e h . 

From straightforward computations, we have the following: 

Av = he h 1 Ari + 2ie h 1 (£ x • Vri) 

d Xj (Av) = he~ r 1 d x . (An) + i(ije ~^An + 2\e~^d Xj (( 2 • Vn) 

2 ix-Cj 

- j^Cije~(Ci ■ Vn). 


ix-<^2 ia3-^2 ix-^2 

d Xk d Xj (Au 2 ) = he~d Xk d Xj (Ar 2 ) + K 2k e~ d Xj (Ar 2 ) + i( 2j e~ d Xk (Ar 2 ) 

1 ia>C2 ia?-C2 

- ■j2 2 j^ 2ke ~ fr Ar 2 + 2\e~h~d Xk d Xj (( 2 ■ Vr 2 ) 

2 ia; -C2 2 i^-C2 

-^<0fce h d Xj (( 2 -Vr 2 ) - -( 2j e h d Xk (( 2 • Vr 2 ) 

-p( 2j (2 t e^((2-Vr 2 ). 
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ix-^2 m 

^xm^ocidxv.dxjU 2 = he h d x dxi^xudxA 7*2 H - 1 C 2 m£ h dxidxudxjT 2 


Tm^Xi^Xk^Xj 1 

^Xk^Xj 1 


ix-^2 1 ice -£2 

+ K2ie~^~d Xm d Xk d x .r 2 - -C,2iC,2me~>^d Xk d Xj r 2 


'Xl^Xk^Xj 1 
^Xk^Xj 1 


ix-^2 icc-C2 

~j^C2mC2k^ h d X id X jT 2 1 C 2 h 9x rn ^xi^Xj'^2 


1 icc-C2 1 icc-C2 

-j~C2i(2ke h d Xrn d Xj r2 j^&m&i&ke h d x ^2 

1 icc-C2 . ia;-C2 

-j^(, 2 mC 2 j& h 9 Xl d Xk r 2 T 1 C 2 h d Xrn d xl d Xk r 2 

1 ix-C,2 1 ia;-C2 

- -j< 2 iC, 2 je~d Xm d Xk r 2 - j^C, 2 m C, 2 iC, 2 je~d Xk r 2 

i ia;-C2 1 icc-C2 

- J^( 2 m( 2 j( 2 ke— d xi r 2 - -C2jC2fc e_fr d Xm d xl r 2 

i icc-^2 1 ia>C2 

- -^< 2 i{ 2 j( 2 k e—d Xm r 2 + ■j^C 2 m( 2 iC 2 jC 2 ke ~ f ^~(1 + hr 2 ). 

Now using the above derivatives, it is straightforward to show the following: 
, ... .. Ce~h ~.... C 2 r 

( 14 ) WMlmm < -^2-lkill^ (n) < ^ h - 

(15) 

(16) 


... .. C 2 R 

||Am 2 ||h2(Q) ^J eh - 

II II / C 2R 

11^2 ||// 4 (Q) < J£ e h ■ 


Therefore we have 


[ d u (A(ui - u 2 ))vdS + f d u (ui — u 2 )(Av)dS\ 

Jan Jan 


C 2R C Mlx,C 2 R C 2 R 

+ — e h )(—e h + —e h j 


-e ft 


<c||v„ JV"«U h4 - ■ A - 

< £e“£e“lW„ -VJ < -M n \ 


1 — 

Now using the fact that ^ < e~h , we obtain 


| f d u (A(ui - u 2 ))vdS + f d u (ui — u 2 )(Av)dS\ < Ce 9 h \\J\f qi - J\f qo \\. 
Jan Jan 

Extending q\ , q 2 to M n by 0 and using (fl2l) and (fT3l) , we get the estimate 


\{ qi -q2m\<C{e-\\M qi -M q2 \\+h) 


Now 


112 II 112 

I'll ~ 1211^-1(0) < ||<?1 — 92||^--1 (r«) 


/ 

J\t\<P 


\(Qi~Q2mw, , 

1 _l_ \C\2 + 


I(« 1 -®)( 0 I S 


\t\>p 


d£, 


1 + ld 2 


1 + ICI 2 







CHOUDHURY AND KRISHNAN 


for appropriate p to be chosen later. 
But 

f |(?T-^?2)(0I 2 ^ / 

12 al = ~ 


l(gi-92)(OP 


f l€l>p 


l + l^l 5 


1 


l£l>p l + P 
2 


dt 


< -o 91 - 92 
P 2 


C 

0-^2 


and 


'lfl<P 


Ifo-ftXOI’ r f 
F“ ^ 1 


(e~ || jV 91 — A/" ?2 1| + /i) 2 


1 + iei 2 


l$l<p 


1 + ICI 2 


■d^ 


<C(e^\\M qi -M q2 f + h 2 ) f - 

J\i\<P 1 


dt, 


+ l€l s 


< Cp n (e~f^\\J\f qi — J\f q2 \\ 2 + h 2 ) 


Therefore 


|?i -®|| ff -i ( n) < C P n e h ll-^9i — d\f q2 \\ +Cp' l h + 


C 


Now assume that ||jV gi — J\f q2 1| <5 = e Then we choose p = {^j^|ln HA/"^ — 
2 

N q2 1||} 71+2 • Further let h = n+2 • With this choice of h, we show that h < ho 

p~^r 

and 1 — /i 2 |£| 2 /4 > 0 for |£| < p. The fact that h < ho follows from these 
inequalities: 

20R 

11 A/"ji -Af q2 \\<e h o < 1 
20 R 


In HA/qrj M q2 \\ < 


hi ||A r qi — A/g 2 HI > 


ho 

2 m 

ho 

1 


^ 20i? |ln||A/gi 

=>- h < ho- 

Now we show that 1 — h 2 —j- > 0 for |£| < p. We have that 

p” = { ^ | to || Y ai -V, 1 |||) ^. 

Since > 1 and 7^r| In — A/" 92 ||| > 1, we have that p n > 1. 
Hence 

*»!«!! <*£ = _L <, 

4 4 4 p n 

and so 1 — h 2 ^- > 0. 
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Therefore 


119 C 18-R n - r 119 4 

||?l-92||jy-i(n) ^ h -^,11 + Chn+2 
h n + 2 

< Ce^\\M qi ~M q2 \\ 2 + Ch^ 2 . 

and since % = t^| In ||jV gi — J\f q2 1| |, we then obtain the estimate 

Iki _ < 72 lll/-i(o) — C{\\N qi — M q2 \\ + | In \\Af qi - A/’qjlH n + 2 ), 

when ||A f qi - M q2 \\ <5 = e“^. 

The case when ||jV gi — N q . 2 1| > 6 follows from the continuous inclusions 

L°°(n) -A L 2 (Q) -a H~\n). 

In other words, we have 


n |, 2 „2 4CM 2 ACM 2 r 

Ikl ^72lliY-i(O) — ^ll^l 92|li,oo(Q) < ^ ^ ||AT gi AU 

and hence the desired estimate follows. This concludes the proof. □ 


5. Stability estimate for slightly more than half data 

Here we prove stability estimates for the partial data case. In the ap¬ 
pendix, we include a proof of the identifiability in this case using linear 
Carleman weights. We would be using a few estimates derived therein in 
this section. 


Proof of Theorem 111. HI We begin with the following identity as at the begin¬ 
ning of Theorem 12.11 and rewrite it as 


(<?2 — qi) U 2 vdx = J d u (A(ui — U 2 ))v dS + j d u {u\ — U 2 ){Av) dS 
i an 

j d u (A(m - u 2 ))v dS + J d v {v,\ - u 2 )(Av) dS 

a^-,e dfl-,e 

(18) + J d u (A(ui - u 2 ))v dS + J d u (ui - u 2 )(Av) dS. 


n 

(17) 


an 


an+,« 


an+,< 
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We estimate the terms in ED- Proceeding as with the full data case, we 
have 

| ( d u (A(ui — u 2 ))v d S + f d u (u\ — u 2 )(Av) dS 1 ] 

Jdft- e Jdri- e 


< C - «2))|| Hi(an _ ie} + R(«l - “2)ll fr 5 (an _ )e) 

(IMUpn) + 

= C\\(AT qi ~ N q2 ){f,g)\\ ^(IMI# 1 ^) + ll^lli/pn)) 

< C\\Af qi ~ Af q2 \\\\(f,g)\\ H i,3^ a ^(\\v\\ H i^ + ||Au||tfi(Q)) 

< C\\Kf qi - ATq 2 \\(\\u 2 \\ H 4(n) + ||A'U2||H 2 (r2))(lklk 1 (n) + 11^^11^1(0)) 

9 R ■—- 

< Ce-^\\M qi -M q2 II. 

Now we estimate the terms in (1181) . 

We first have 

| f d u (A(u\ — u 2 ))v d5| = | f e ^d v {A{u\ — u 2 ))e~v AS\ 

Jdfl- i_ e Jd£l+ e 


+ ,e 

x ■ a. 


< lie h d v (A(ui - u 2 ))\\ L 2 {dn+ f) \\e * u|| L 2( 90+ e) 

< C\\e ~d u (A(ui - u 2 ))\\ L 2 (dn +tt ) 


and 


L 


d u (u\ — u 2 )(Av) d5| = | / e h d u {u\ — u 2 )e h (Av) dS| 


< ||e h d u (ui — u 2 )\\ L 2( dn+ ^\\e h (Av)\\ L 2^ dQ+ ^ 

< C\\e ~d u (u\ - u 2 )\\ L 2 ( d n + e y 
By the boundary Carleman estimate, we have for e > 0, 

h*\\y/a ■ v eT TT <9„(-/i 2 Au)|| L 2 (an+) < C\\\e~^{h A B qi )u\\ L 2 {Q) + 
h^\W~a ■ v e~^d v (-h 2 Au)|| L 2 (90 _) + /i§|| y/-{a ■ v) e _ ^^u|| £ ,2 (9n _ ) ) 

We then have 

II y/a • v e~ 2 ^d u (Au)\\ L 2 {dn+) <c(Vh\\e~ li ^B qi u\\ L 2 ^ ) 

+ ^\\V~a ■ v e _ TT^u|| i 2 (9n _ ) 

+ || \J-ol-v e~^d u {Au)\\ L 2 m _A. 
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Since on d£l +iS , a ■ v > e, we have 


e x h d v {Au)\\ L 2 {m+e) <^=(Vh\\e B qi u\\ L 2 {n) + \\W-{a ■ v) e V^u|| L 2 (9n _ 


+ II y/-{a ■ v) e“^^(A'u)|| L 2 (an _ ) ). 

< ^(yh\\e~^B qi u\\ L 2 m 

+ J- jnf (a ■ iy)\\e~^d u (Au)\\ L 2 {dn _ E) 


+ J-jnf (a-p)||e h d u u\\ L 2 ^ m _^j. 


Therefore, 


\[ d v (A(ui - u 2 ))v dS| < -^=, /- inf (a ■ v) \\e *h d v (/\(ui - u 2 ))\\ L 2 m _ \ 

Jdn+, e Vey 


+ C\/-\\e h (B qi )(ui-u 2 )\\ L 2 (n) 


C 


+ Jh/ey sn f ( a ' ^ h ^( Ul-M2 )H R 2 (an_, £ )- 


Now we have 


||e x h B qi (ui-u 2 )\\ L 2 {n) = \\e V ( gi - g 2 )u 2 || L 2 (n) . 
Using the CGO solutions from Proposition (13.31) . 


i^<2 


u 2 (x, ( 2 ; h) = e h (1 + hr 2 (x, C 2 ; h)) where 




lei 5 


C 2 = —7T + \/ 1 - b 2 —rP - i 


ice, 


we have 


|e X h( qi - q 2 )u 2 ||z, 2 (n) < C. 


Therefore 


ian +] 


d„(A(ui -u 2 ))v dS| < CVh + Ceh\\d u (A(ui - u 2 ))\\ L 2 ^ 9n _^ 


C R 

+ -^ e ~ K \\dv{ui -u 2 ) ||l 2 ( 9 o_, £ )- 


R 

eh 


<C[yh + —( ||ft,(A(«i - «2))|| £ 2( an _ ie ) 
+ \\d u {ui — «2)||L a (9n_, e )))- 


C 


< CVh+-j~ eh II Nqi - A/’g 2 ||(||u2||_H 4 (n) + l|Au2||lL 2 (Q)) 
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and using the estimates (p~6li and it follows that this is 


< CVh + —^e h || J\f qi — M q2 \\ 

< CVh + Ce h \\Af qi — A fq 2 1|, 

where the constant C now depends upon e. 

From the boundary Carleman estimate, we also have 

h?\Wa ■ v e~^d u u\\ L 2 idn+ e) <C(\\e~^(h 4 B qi )u\\ L 2 ^ n) 

+ h*\\yj-a-v e~ 2L ^d u (-h 2 Au)\\ L 2 ^ gn _ t£) 

+ h%\\y/-a-u e-^d u ii\\ L 2( dn _ t£) ) 

and therefore 

IIa/ a ■ » e~^d u u\\ L 2 (gn+ e) <c[}i^\\e-^B qi u\\ L 2 {n) 

+ h\\y/-a ■ v e~^d u (Au)\\ L 2 {gn _ !£) 

+ \\V~( a ■ v) e~^d u u\\ L 2 {gn _ y 


We then have 


C 


\e h d v u\\ L 2 , gn ) < —f(/i 2 ||e h B qi u\\ L 2 , n) 

V e 

J rh\\^-(a ■ v) e~^d p (Au)\\ L 2 {gn _ t£) 
+ || V-(a • v) e~^d u u\\ L 2 ^ gn _ j£) ) 

< ^j=h^\\e~^(B qi )u\\ L 2 {n) 
v e 
hC 


h ( I x . a 

+ ~ (“ ' V )\\ e ~ — d ^U)\\ L 2 {aQ _^ 


C 


+ 7ey Mi a ' l/ ll e h d " u W L 2 (dn.,e)- 
Similar to the previous estimate, we have 

f _ C 3 g-a 

| / d v (ui - u 2 )(Av) dS\ < —=h^ ||e h (B qi )(ui - u 2 )||L 2 (f 2 ) 

Jda+, t \/e 

hC 1 i - 

+ ^l-Ma-v\\e-^d u (A( Ul - u 2 ))|| La(fln _ >1 
C I - 

+ ~r\h io f ( a ' ^ W e ~~^M u i “ «2)||L2(an_, E ) 
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where the constant C again depends upon e. 

Therefore using the estimates obtained above, we have 

| [ e~ lx '^(q 2 — <71 ) dx\ < \ f d u (A(ui — U 2 ))v dS\ + \ f d v {u\ - U 2 ){Av) dS\ 

J ft J d£l-J 

+ \[ d u (A(ui - U 2 ))v dS + [ d v (u\ - u 2 ){Av) dS| 

+ 1 / (<?2 - qi)e~ lx ^(hri + hr 2 + h 2 fir 2 ) dx\ 

Jn 

< CVh + Ce™ \\Kf qi -Kf q2 \\+Chi 

+ Ce h ||jV 9l — M q2 1| 

9 R 

+ Ce h ||A f qi — A/^ 2 1| + Ch 
( 19 ) <C(Vh + e^\\Sf qi -Sfq 2 \\). 


The argument that now follows is similar to the one in [Sj. We will apply 
Vessella’s result given in Theorem 13.41 for the following set up. We take D 
to be the ball 5(0, 2) and E = V PI 5(0,1) where V is a suitable small open 
cone centered at 0 obtained by perturbing the vector a. slightly and recalling 
that £ is perpendicular to a. Note that the above estimate is valid for all 
£ € V such that |£| < -|. 

Now let q = q\ — q 2 extended to M n as 0 outside and for a fixed p E (0, |), 
let /(£) = q(pO- Then / is analytic in 5(0,2) and 


P np 

T”/«)l < s 2M l n l° l ( d iM .(«)-)i 


-iOal ' 


Taking C and A in Vessella’s result to be C = 2M\Q\e np and A = diam(fl) 1 , 
we get that there exists a constant 71 € ( 0 , 1 ) such that 


|/(£)| < (7 1 -T' 1 ( \E\/\ D \) 


sup 1/(01 

E 


7i(|B|/P|) 


for all £ € 5(0,1). 


Letting 0 = 71 1571/1ZA|, we have that for all |£| < p, 


( 20 ) 


15(01 < c 


l-e 


sup |g(£/p)| 

K vnB(o, P ) 


Note that the constant 8 is independent of p and h. We have 


q 


2 

e 




([ I5(0I 2 

Vw + ki 2 


d£ + 


r \m \ 2 

h\> P 1 + l^l 2 



1 

e 


< c 



2 

e 


L°°(B(0,p)) 
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The estimate of the second term on the right hand side above is obtained 
from Plancherel identity. Now from f)20[i . it follows that the left hand side 
is 


<C^pee 2np o lk(OllL°°(RnB( 0 ,p)) + 2 


<C(^e 2 ^ e f \W qi -N, 


'92 


I 2 + p n e e 2npl 8° h + \ | . 


1 

Let L = 3n + 2 ~ 2n9 and 5 = e ~ eK/h ° , where K = ^±2 + 4n 1=2 + i 8 R. 


Let \\AT qi - AT q2 \\ < S. 

Then choose p = In | In \\M q ^ — J\f q2 1| | and 


h = 


n+2 2np(l — 9) 

p o e o 


Claim 1: p < j . We have ph = 


1 


n + 2 2np(l —0) 

)S~e - 0 - 


< 


1 


n+2 2np(l — 0) — n + 2 • 

p~S~ 1 e 0 P^^ 

Now since n > 3 we have — 1 > 4, and hence p~ 1 > p 4 . Therefore, 
ph < jp < 2 (since p > 1). 

Claim 2: h < ho- We have 


||A/"gi M q2 1| < e 


i 

-e K ' h 0 


< 1 


In H-A/^r, -N q2 1| < -e K//l o 


In ||2V Q1 - 1| | > 


K 


In | In ||jV 91 -M q2 HI > 


hi 


=> h = 


^ p^ > 

3n+2—2n6 


> 


ho 

1 


1 


< 


n+2 2n(l-9) 

P 8 p 9 

1 


ho 

< ho 


n+2 2np(l — Q) — n+2 2n(l — 6) 

p o e o pope 


< K 


o 


Claim 3: 1— h 2 ^- > 0. This is because h 2 — j- < h 2 - j- = 


2 <" + 2 > 4nptl^i 


2 — 2 yi +4 


4p 0 e' 


„2— 


2n+4 


4e' 


,4np 


f=g- < p* s = 2 n+ 4 —• Now 2Ti n t ~ — 2 > 8 since n > 3 and therefore 
























STABILITY ESTIMATES FOR A BIHARMONIC BOUNDARY VALUE PROBLEM 15 


2n+4 _o o - . - . 1 

p e > p which m turn implies that 


h 2 


i£i 2 i 

< -=■ < i. 


Then we have 

f + = n ^[2npi^+18 R(p 


n-\- 2 


p & e 


h = p6 e l- 


2np- 


i a n+2 o_ 1 — 6 

< e fP+2npV+i8«P^- 2rip ^ since p* < 

< g (f P+2npL r ^)_|_ e ( 1 8H+2^p+2npL ;r ^) 

o (f p+2npi ¥ ^) (i8iip+ll^p+2npT g A) 


< e 


since p > 1 


/ n , r, 1 — 0 i i o r> i n+2 , r, 1 — 0 \ 

( -n-+2n— 7 j |- 18 i 2 H 75 K2n—75— )p , . . 

< Ce eK ’ since e a + e b < 1 + e a+b . 

Therefore since K = 2w g f2 + 4n^p + 18 R and p = In | In ||A/" gi — jV 92 |||, 
we obtain 

Iklln-ifni - C ( 11-^91 -^ 92 II + (-^bilM-A/gi — A^ 2 lll) - ^) 


and hence 


liL-i(H) < C(\\AT qi - Af q2 \\ + (— In | In \\J\f qi - J\f q 2 \W) 0 ) 2 , 


whenever ||A f qi — M q2 \\ < 5. 

When \\J\f qi — J\f q2 \| >5, we have 

„ 2 CM e 2CM ~ ~ M e 

II?i — Q2\\H- 1 (n) < C 7 ||gi — <?2||L°°(f2) < — —b 2 < —— 1 | M qi — -Wj 2 || 2 

62 5 2 

and the desired estimate follows. 


□ 


Appendix A. 

In this section we prove the unique determination of q from (J2J) when the 
Neumann data N q is known on slightly more than half the boundary. This 
is already done in a more general set-up with limiting Carleman weights in 
eh , where the authors use logarithmic weights. We give here the proof with 
linear Carleman weight following [2\ for the sake of completeness. 

Theorem A.l. [11] Let II C R n ,n > 3 be a bounded domain with smooth 
boundary. Consider Equation © for two potentials qi,q 2 £ L°°(Q). Let 
J\f qi and J\f q2 be the corresponding Dirichlet-to-Neumann maps measured on 
If N qi = N q2; then qi = q 2 . 

Proof. As before, we start with the following integral identity. 

/ (Q 2 —Qi)u 2 V dx = — / d u (—A(ui~u 2 ))v dS— / d v (u\—u 2 ){— Av) dS 

Jfl J 
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L 


d u (A(u\ — U 2 ))v els’! = | / e h d^(A(ui — U 2 ))e h v dS| 
90+,,= Jdfl+, e 


L 


< ||e X hd v (A( Ul -u 2 ))\\ L 2 {dn+!e) \\e X h V \\ L2{m+e) 


d v {u i — U2){Av) d^l = | / d v (u\ — U2)e h ( Av ) dS^ 

0O_)_ ;<E J dfl .|_ j£ 


< ||e h a i/ («i-« 2 )||L 2 (an+ie) ||e h (Av)|| i 2 (sn+ie) 

From the boundary Carleman estimate, we have 

Q 

\W a ' v d v (-h 2 Au)\\ L 2 [dn+) < — \\e-^(h A B qi )u\\ L 2 {n) 

h 2 

where u = u\ — U 2 - 
Using this we get 

V'eWe-^d u {-h 2 Au)\\ L 2 {dn+ e) < ||V« ■ v e~^d v {-h 2 Au)\\ L 2 {dn+ t) 

< ||V« • v e~^d v (-h 2 Au)\\ L 2 ( dn+) 

C 


< —||e h (h B qi )u || L 2 (q) . 

h 2 


Therefore 


|e V^(A(txi-« 2 ))||L2(an + , e ) < c '\/^ll e ^ fe - 9i)e h 2 (l + /ir 2 )|| L 2 (f2) 


< C 


+ ^ r 2)||L 2 (D) < C 


From the boundary Carleman estimate, we also have 

(j 

\\y/ot ■ v e~^d u u\\ L 2 {dn+) < - T ||e” _ Jr(/i 4 B, 1 )i4|| £ 2 (n) 

h 2 

Again using this, we get 

Ve||e~ 2 ^3 1/ iz|| L 2 (an+i£) < || V ol ■ ^ e -a 7rd„ii||2^(oi +> .) 

< || \faTv e~^d u u\\ L 2 {dn+) 

< ^||e _ ^(/i 4 i3 g i)u|| L 2 (a ) 

/i 2 

(7 

=>- ||e _:2 ^^(wi -tt2)||L2 ( an+ e ) < _ ^ — zll e_2 '^(^ 4 ^gi) u IU 2 (n) 

h2y/e 

Ch^ ... , ... Chi 

< ——11(^1 + /ie 2 )|| L 2 (r2) < 






Next we show that the terms ||e J* 'e’|| J c 2 (on_ t _ e ) and ||e (Au)|| i 2(oi+ \ are 

bounded. The term 

l|e-^w||L2(an +ie) < ||e"K"i;||i2( a n) = ||(1 + hrT)||L 2 (dn) 

<C(l + ||hr 1 || Hl(n) )<C, 
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since li<l. 

Again 

(Ax) = h ArT- 2ie ^ ( “ ^'0 ^ Vn. 

Therefore, we have 


||e~(A-u)|| L 2 ( an+je ) < \\e~ (Ax) ||z, 2 (V 2 ) 

< \\hAri\\ L 2( d ty + C\\Vri\\ L 2( d ty 

< C(h\\Ari\\ H i(ty + ||Vri|| H i(Q)) 

< c {-^\\ri\\ H tJQ) + ^lki||^4 i( n)) < C. 

Also using the estimates on ri,r 2 it follows that as limit h —> 0, 

/ (?2 - qi)u 2 v dx -A / e~ lx< {q 2 - qi) dx. 

Ju Jfi 

Therefore combining all the above estimates and passing to the limit as 
h—>■ 0, we have 


/ 

Jn 


e lx '^{q 2 ~ qi) dx = 0 


for all ^ € M n perpendicular to a. Varying a in a sufficiently small neigh¬ 
borhood, we see that above estimates is true for all £ in an open cone in R n . 
A simple application of the Paley-Wiener theorem then implies that q 2 = q\ 
on 17. This concludes the proof. □ 
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